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Abstract. 

We investigate some simple and surprising properties of a one-dimensional 
Brownian trajectory with diffusion coefficient D that starts at the origin and reaches 
X either: (i) at time T or (ii) for the first time at time T. We determine the most 
likely location of the first-passage trajectory from (0,0) to {X,T) and its distribution 
at any intermediate time t < T. A first-passage path typically starts out by being 
repelled from its final location when X 2 /DT -C 1. We also determine the distribution 
of times when the trajectory first crosses and last crosses an arbitrary intermediate 
position x < X. The distribution of first-crossing times may be unimodal or bimodal, 
depending on whether X 2 /DT < 1 or X 2 /DT 1. The form of the first-crossing 
probability in the bimodal regime is qualitatively similar to, but more singular than, 
the well-known arcsine law. 


PACS numbers: 05.40.Jc, 02.50.-r, 05.40.-a 
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While Brownian motion has been extensively studied and much is well understood about 
this process (see, e.g., M), a number of simple questions continue to be investigated 
when the motion is subject to a constraint. Some examples of this genre include 
determining the time that Brownian path spends in a range dx about a point x the 
area swept out by a Brownian motion when it first returns to the origin ms, the time 
when a one-dimensional Brownian motion, which starts at x > 0, attains its maximum 
before first crossing the origin [T2j[T3], as well as related extremal properties of such 
motions m- There has also been considerable study on the topic of level-crossings, 
namely, the crossings of a given point on the line, for both Brownian motion [ 151116] 
and for general stochastic processes mm, as well as investigations of the crossings of 
a time-dependent level (see, e.g., 


(X,T) 



Figure 1. Illustration of the quantities under study. Shown is a schematic first-passage 
path from (0,0) to (X,T). At some intermediate time (dashed line), the particle is at 
position x. At some intermediate level (dotted line), the first- and last-crossing times 
of this level are tf and te. 

Given these extensive studies, it is surprising that some basic issues about level 
crossings seem to have not yet been fully addressed. In this work, we investigate the 
intermediate-level crossings in both space and time of a Brownian motion that first 
reaches a final destination X > 0 at time T. This type of question underlies the task 
of reconstructing a Brownian path from measurements at discrete times, an issue that 
is of basic importance for understanding diffusion in porous media (see, e.g., [25H2Z]). 
The condition that the particle must reach X for the first time at time T constrains the 
entire trajectory before time T in a non-trivial way. We will study the consequences of 
this constraint on the following basic properties of a Brownian particle with diffusivity 
D at intermediate stages on its way to its first passage at (X, T), as illustrated in Fig. [H 

(i) At an arbitrary intermediate time, what is the expected position x of the particle? 

(ii) When does the particle first cross an intermediate level x < XI 

(iii) When does the particle last cross this intermediate level? 
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Our approach to answer these questions is based on decomposing the first-passage 
probability to (X, T ) as the convolution of the propagator to the intermediate point 
(x,t) and the propagator from this intermediate to the final point. In this framework, 
the intermediate crossing properties are simple to formulate; however, the results are 
somewhat surprising and depend in an important way on the global nature of the full 
first-passage path. In the “ballistic limit”, where X 2 3> DT , the space-time trajectory 
of the first-passage path approaches a straight line. In this case, intermediate-crossing 
properties are close to those that arise by treating the first-passage trajectory as purely 
ballistic. In the opposite “diffusive limit”, where X 2 <C DT , the first-passage path is 
typically repelled from its final location X so as to avoid hitting X before time T. 

Correspondingly, the probability distribution of times for the particle to first cross 
an intermediate level 0 < x < X undergoes a transition from unimodality, for X 2 DT , 
to bimodality, for X 2 -C DT. In the latter case, the first-passage path to ( X , T) is most 
likely to cross an intermediate level near the beginning or near the end of its trajectory. 
In this bimodal regime, the distribution of the first-crossing times is reminiscent of, but 
more singular than, the famous arcsine law [3j[3] for the distribution of times that a one- 
dimension Brownian particle spends on one side of the origin, as well as the distribution 
of the last passage time at the origin. 

As a preliminary, we first investigate intermediate crossing phenomena for a purely 
Brownian path between (0, 0) and (X, T) in section [2] In spite of the simplicity of 
this problem, the distributions of the first- and last-crossing times of an intermediate 
0 < x < X level are also quite rich and have a qualitatively similar unimodal to 
bimodal transition as in the case where the trajectory is a first-passage path. Moreover, 
the boundary in phase space that demarcates this transition is, surprisingly, not single 
valued. In section [3j we investigate intermediate crossing properties when the trajectory 
between (0, 0) and (X, T) is a first-passage path. We briefly discuss the situation where 
the intermediate level is negative in section |4j Finally, we offer some perspectives in 
section [5j 

2. Preliminaries: Intermediate Crossings for Unconstrained Brownian 
Motion 

Consider a Brownian particle in one dimension that starts at the origin. We will 
make extensive use of two important characteristics of this Brownian motion: (i) 
the occupation probability and (ii) the first-passage probability. These quantities are, 
respectively: 

™ = vmr e ' XVlDT F(A - r) = A e ' XVtDT ' (1) 

Here P(X, T) is the probability that a Brownian particle moves a distance X from its 
starting point over a time T, while F(X, T) is the probability that the particle first 
reaches this point at time T |29j . One of the intriguing aspect of a Brownian trajectory 
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in one dimension is that it ultimately reaches any point on the infinite line. However, 
even though the trajectory is guaranteed to eventually reach any point X ^ 0, the 
mean time to reach an arbitrary point is infinite. This dichotomy arises because the 
probability that the trajectory first reaches position X at time T asymptotically decays 
as T~ 3 / 2 , from which the mean time to reach X is infinite. 

For a Brownian particle that starts at (0, 0) and ends at ( X , T), we now ask: where 
is the particle and what is its probability distribution of positions at an earlier time 
t < T? The probability that the particle propagates from (0,0) to (X, T) may be 
written as the convolution of the probabilities for the path from (0, 0) to ( x , t ) and the 
path from ( x,t ) to the final location (X, T): 


P ( X,T ) -fp^)PiX- X ,T-t )dX 

r p -x 2 /ADt p -(X-x) 2 /AD(T-t) 

= / 7 _ -^=^=- dx. (2) 

J yfteFDt ^4vr D(T-t) 

Integrating over x reproduces, after some simple algebra, P(X,t ) in Eq. (JTJ). For later 
convenience we introduce the dimensionless time r = t/T and the dimensionless length 
X = x/X. By maximizing the integrand with respect to x. it is straightforward to derive 
that the integrand has its maximum when x = tX. That is, the particle moves a 
fraction r of the final distance in a fraction r of the total time, a well-known property 
of Brownian motion [281] . 

Finally, we compute the probability V(x,t) that the particle is at x at time t, 
given that it is at X at time T. This quantity equals the probability for the subset of 
Brownian paths from (0, 0) to (X, T) that reach x at time t divided by the probability 
for all Brownian paths from (0,0) to (X, T): 


V{x,t ) 


P(x,t) P(X—x, T—t) 

nxTr) 

£ -x 2 /ADt e ~{X~x) 2 /AD{T-t) 

y/AnDt ^AnD(T-t) 

e -(x-rX) 2 /ADTt(1-t) 


x VWPe x2 / 4DT 


4:TtDTt(1 — t ) 

e -a(x-H 2 A(l-r) 

7T vM 1 “ T ) 


(3) 


where we have rewritten the last line in terms of y e |!(—oo, 1), r e (0,1), and the 
dimensionless parameter a = X 2 /ADT , which demarcates whether the full trajectory is 
ballistic, for a 1, or diffusive, for a 1. 

From Eq. (|3jl. the probability distribution to reach an intermediate point x at 
dimensionless time r is a Gaussian centered at y = r whose width \J r(l — r) is maximal 
for t = | and vanishes for r —> 0 and r —> 1. The physical message from these results is 
that the average position of a Brownian particle that starts at (0, 0) and ends at (X, T ) 
interpolates linearly between 0 and X as the time increases from 0 to T [28] , 
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In a similar vein, we now determine when the Brownian particle first crosses and last 
crosses an intermediate position x < X, given that the particle is at X (not necessarily 
for the first time) at time T. More generally, we compute the probability distributions 
for these two events. Using the same reasoning that led to Eq. (J3]), the probability 
T(x,t) that the particle first crosses x at time t is given by: 


T(x, t ) 


F(x, t) P(X — x,T — t) 
P(X,T) 


(4a) 


That is, the first-crossing probability equals the probability for a Brownian particle, 
which starts at (0,0), to first reach (x,t) —hence the factor F(x,t )—times the 
probability that the particle propagates from (x,t) to (X, T), normalized by the 
probability for all Brownian paths that propagate from (0,0) to (X, T). 

Similarly, the probability jC(x,t) that the last crossing of x occurs at time t is: 


C(x, t ) 


P(x, t) F(X — x,T — t) 
P(X,T) 


(4b) 


The reasoning that underlies (14bjl is slightly more involved than that for (Uajl . For (x, t ) 
to be the last crossing, the remaining trajectory must be a time-reversed first-passage 
path from (X, T) to (x,t). This constraint guarantees that the crossing at x is the last 
one in the full trajectory to (X, T). 

Substituting the explicit forms for the occupation and first-passage probabilities 
from Eqs. (OP) and performing some straightforward algebra, the first-crossing and last¬ 
crossing probabilities are, after expressing all variables in dimensionless form: 


P(x,r) 

£(x , T ) 



g —a(x—r) 2 /[r( 1 _ r)] 
e -a(x-r) 2 /[r(l-r)] 


(5) 


It is apparent that X(y,r) = £(1 — y, 1— r), and vice versa, so that we only need 
to study one of these quantities. Figure [2] shows the time dependence of the first- 
crossing probability P(x, T ) for three representative values of a and y. The first-crossing 
probability clearly becomes bimodal for small a and y —y 1. On the other hand, for 
a —> oo, so that the trajectory becomes progressively more ballistic, the Erst- and 
last-crossing probabilities must approach each other. 

To determine where in the y-a plane that the first-crossing probability changes 
from unimodal to bimodal, we first need to locate the extrema of J 7 ( y, r) for each value 
of y. Differentiating X(y,r) with respect to r and setting the result to zero gives the 
cubic equation 


4t 3 + (2a — 4ay — 7)r 2 + (3 + 4ay 2 )r — 2ay 2 = 0 , (6a) 


and solving for a gives 


(4 t 3 — 7 t 2 + 3 r) 

[(4y - 2)t 2 - 4y 2 r + 2y 2 ] ' 


a = 


(6b) 
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X = 0.2 X = 0.5 x = 0.8 



0.0 0.5 1.0 0.5 1.0 0.5 1.0 


T 


Figure 2. First-crossing probability _F(x,t) (red, solid) and last-crossing probability 
£(X,t) (blue, dashed) for illustrative values of a = X 2 /ADT and x = x/X when the 
initial trajectory is an otherwise unrestricted Brownian path from (0,0) to (X, T ). For 
small a, F{X: T ) changes from unimodal to bimodal as x increases. For large a (almost 
ballistic trajectory), the first-crossing and last-crossing times are close to each other. 




(a) 


(b) 


Figure 3. (a) Loci of the extrema of .F(x,t) in the a-r plane for various values of 
X when the initial trajectory is an otherwise unrestricted Brownian path to (X,T). 
For a given a and x, the extrema occur at the values of r where the line a = const, 
intersects the curve that corresponds to a given x- (b) Zoom near the critical point 

(f,l) 
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These loci for a are plotted as a function of r for various values of x in Fig- [Ufa). For 
given a and x, the extrema occur at the r values where the line a = const, intersects 
the branches of the curve for a given x- For example, in Fig. [3](a), when a = 0.5 
and x — 0.9, T has three extrema at r = 0.32, 0.64, and 0.99 (two maxima and an 
intermediate minimum), as indicated by the dots. On the other hand for a = 1.25 and 
X = 0.7, there is a single extremum (circle). When x ~ 0.75, the regime of bimodality 
extends over the widest range of a. We plot the phase boundary (Fig. [jj on the x~ a 
plane by scanning across x and numerically finding the range of a values at which three 
extrema occur. 




(a) 


(b) 


Figure 4. (a) Numerically obtained phase diagram from Eq. ([6]), showing the domains 
in the parameter space where the first-crossing probability is unimodal and where 
it is bimodal for an unrestricted Brownian path from (0,0) to (X,T). The points 
marked a,6,c are at (0,| — 2-\/3), (0.747835-•• ,1.09325...), and (1, |), respectively, 
(b) Detail near the cusp in the phase boundary; the points marked d, e are at (|,1) 
and (|,4(2 —-\/3)), respectively. 


The behavior near the critical point (r, a) = (|, 1) is particularly intriguing 
(Fig. |3](b) ). For a given value of x that is less than but very close to x = 3/4, there are 
two distinct sets of solutions—one set with a slightly less than 1 and another set with 
a slightly larger than 1. There is also a small gap between these two solution sets, as 
indicated in Fig. [3](b) for the cases x = 0.749 and x = 0.7499. The consequence of this 
feature is that the unimodal to bimodal phase boundary is not single valued near the 
cusp, as shown in Fig. |4j It is quite remarkable that such a simple question about a 
Brownian path—namely the time for the path to cross an intermediate position—leads 
to such a rich phenomenon. 


3. Intermediate Crossings for First-Passage Paths 

We now turn to the properties of intermediate crossings for first-passage paths, where 
the first-passage constraint for the full trajectory from (0,0) to (X, T) significantly 
affects the global nature of intermediate trajectory and concomitantly, the intermediate 
crossings. We divide our discussion into: (i) the location of the particle at an arbitrary 
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intermediate time and (ii) the time when the particle crosses an arbitrary intermediate 
position. 


3.1. Intermediate Crossing Position 

At a time t < T, we again ask: what is the probability V(x,t) that the particle 
is at position x, given that the particle first reaches X at time T? We determine 
this probability, following the same approach that led to Eq. (EJ), by decomposing the 
trajectory into the segment up to time t and the remaining segment from T — t to 
T. The first segment is a Brownian path to (x,t), but now subject to the constraint 
that it can not reach X, because reaching X would mean that the full trajectory is not 
a first-passage path to (X,T). To implement this constraint we impose an absorbing 
boundary condition at X by augmenting the initial Gaussian with an anti-Gaussian 
image distribution that is centered at 2X [31129] . We denote the probability distribution 
in the presence of this absorbing boundary condition as Pa{x,P)- The second segment 
from T — t to time T is simply a first-passage path between these two points. 

Assembling these elements, we have 


V{x ,t) = 


Pa{xC) F(X—x,T—t) 


F(X,T) 


1 


y/AxDt - 


e ~x 2 /4Dt _ ^—(2X—x) 2 /4Dt 


\X — x\ 


d -(X-^c) 2 /4D(T-A) 


j ^4:TTD(T-t) 3 


\X\ 


VAnDT 3 


j—X 2 /4DT 


(7) 


By expressing all quantities in terms of the dimensionless variables introduced 
previously, the above expression simplifies to 


V{ X ,r) 


[a 1 1 ~ XI 

V ^ \/t{ 1 - t) 3 


g-“ [(!—x) 2 —(1—-r)] /(1 -t) 


e ~otx 2 /T _ e -OL(2-x) 2 /T 


f® 11 ^ I e -<*(X-T ) 2 N _ g-4a(l-x)Al 

V 7T \Jr(l - r) 3 


( 8 ) 


Notice that r lies in (0,1) by construction, but \ can range from — oo to 1. 

The behavior of V{Xi T ) i s quite rich, as illustrated in Fig. [5] For small a (diffusive 
limit), the particle must initially move in the negative -x direction so as to ensure that it 
does not hit X before time T. A related type of repulsion arises in the positions of a set of 
N random walkers that all start a finite distance away from an absorbing point subject to 
the constraint that none of them are absorbed up to a given time [30]. As a is increased, 
corresponding to a more ballistic and hence more deterministic trajectory from (0, 0) to 
(X, T), the spread in the probability distribution correspondingly decreases. 

The extent of the repulsion may be quantified by the minimum value of x mp (r), 
the most probable location of the trajectory as a function of rescaled time. For a —> 0, 
the difference between the initial and final positions becomes immaterial so that we can 
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Figure 5. Heatmaps of the conditional occupation probability P(x,r) for a first- 
passage path from (0,0) to {X,T) for illustrative values of a in the r-y plane. 
The smooth curve shows the time dependence of the most probable location of the 
particle. For small a, the trajectory is initially repelled from its final destination. The 
distribution is most diffuse near r = i and most localized near r = 0 and r = 1. 

approximate the minimum value of x mp (r) by x mp (r — \). Using this simplification in 
Eq. (0, the probability distribution at r = | reduces to 

V(y, \) oc y/a e~ a y e~ Aay2 sinh(4m/), (9) 

where y = X~ 1- Finding the maximum of this latter expression is now elementary and 
we find that for a —> 0, the minimum on the most probable trajectory is located at 

U* - ~iTr ■ ( 10 ) 

2y a 

Thus for a —> 0, the most probable trajectory is strongly repelled from the final point 
for r < 1/2 and subsequently is strongly attracted to this final point. We note that this 
behavior for y* can also be inferred from a macroscopic fluctuation theory, as given in 
Ref. [30] . 

3.2. Intermediate Crossing Times 

Parallel to the discussion given for the pure Brownian path, we now compute the 
distribution of times when the first-passage trajectory to ( X , T ) crosses an intermediate 
level x G (0,AT). We shall consider the case x < 0 separately in section [4] as the 
calculation details are different. We also separately consider the cases of the first crossing 
and the last crossing because the calculational details and the results for these two cases 
are quite different. 

3.2.1. First-crossing probability As in the case where the full trajectory is an 
unconstrained Brownian path that starts at (0,0) and ends at (X, T), we decompose 
the trajectory into an initial segment that starts at (0, 0) and crosses x for the first time 

















Intermediate-Level Crossings of a First-Passage Path 


10 


at time t, and a final segment that starts at (x, t ) and crosses X for the first time at 
time T. We consider here the case x > 0; the case x < 0 will be treated in the next 
section. Here the analog of Eq. (Hall is 


F(x, t ) 


F(x, t) F(X — x,T — t) 

FVCT) 


( 11 ) 


We now substitute the relevant first-passage probabilities from © into the above 
equation to obtain 


F(x,t) = 


\x\ 


y/ An Dt 3 


0 -x 2 /ADt_ 


\X 


X 


y/AnD(T-t) 3 


-(X-xf/A D(T-t) 


\X\ 


y/AnDT 3 


d -X 2 /ADT 


( 12 ) 


In terms of the dimensionless variables a, y, and T i the above expression simplifies, after 
some straightforward algebra, to 


F{x,t) 


|1 ~ X\ IXI e -a(r- X ) 2 /h(l-r)] 

V * [( l _ r ) r ] 3/2 


(13) 


Once again, the qualitative behavior of this first-crossing probability depends in an 
essential way on the value of a (Fig. [6]). For large a (ballistic limit) T{x,t) is sharply 
peaked about the point r = y. Thus in this ballistic limit, the most probable time 
for a Brownian particle to first reach a distance x is simply equal to x. Moreover, the 
distribution in (fT3|) reduces to a Gaussian form for r ~ y. 

Conversely, in the diffusive limit of ck —y 0, and for r a and 1 — r 3> a 
(i.e., r not too close to 0 or 1), the probability density function is controlled by the 
factor [r(l — r)] 3 ^ 2 , which is similar to, but more singular than the arcsine law [4]. 
The surprising result for this limit is that a Brownian particle is most likely to cross 
an arbitrary intermediate level either at the very beginning or at the very end of its 
trajectory. 

A remarkable aspect of the first-crossing probability (1131) is its invariance under the 
simultaneous interchanges r —> 1— r and x 1 — X- We can give a simple graphical 
argument to justify this symmetry (Fig. 0). In (a), a first-passage path from (0,0) 
to (X, T ) is comprised of a first crossing (dashed) to (y, r) (in scaled units) and the 
remaining segment to (X,T) (solid). Interchanging these two segments leads to a first- 
crossing segment to (1 — y, 1 — r) and the remaining segment to (X,T). Since the 
segments are independent, the probability for these two first-passage paths in the figure 
are identical and thus X(y,r) = J r (l —y, 1—r). 

We can apply this same perspective to argue that F is bimodal for a 1. Indeed, 
let us compare a first-passage path to the final point that has its first crossing close to 
r = 0 or r = 1, and a first-passage path that has its first crossing near r = 1/2. In the 
former case (see Fig. 04)), the remainder of the path away from the first crossing must 
be repelled from the final point. The probability for such a first-passage path from (0, 0) 
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Figure 6. First-crossing probability F{Xi T ) (red, solid) and the last-crossing 
probability £(y,r) (blue, dashed) for illustrative values of a and x when the initial 
trajectory is a first-passage path from (0,0) to (X, T). For small a, F(x, r) is bimodal 
while C{Xi T ) is unimodal and sharply peaked as r —> 1. For large a (almost ballistic 
trajectory), the first-crossing and last-crossing times nearly coincide. 

(X,T) (X,T) 




Figure 7. (a) A first-passage path that consists of a solid and dashed segment with a 
first-crossing at (y, r). Interchanging the order of these segments gives a first-crossing 
at (1 —X, 1 —t). 

to ( X, T ) scales as T~ 3//2 for a <C 1. On the other hand, the probability for a path that 
first-crosses x near T /2 is the product of the probabilities of two first-passage paths of 
duration T/2, namely (T” 3//2 /2)“, which is much less than T~ 3//2 . Thus a first-crossing 
near T/ 2 is unlikely for a -C 1. 

To determine where in the phase space the transition in T between unimodality and 
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Figure 8. (a) First-passage paths that first cross an intermediate level near r = 0 or 
r = 1 on an exaggerated scale to emphasize the limit a<l. (b) A first-passage path 
that first crosses the intermediate level near r = 1. 


bimodality occurs, we again determine the extrema of T. Following the same procedure 
that led to Eq. (J6]), we obtain the cubic equation for the location of the extrema 

6t 3 + (2a — 4ay — 9)r 2 + (3 + 4ay 2 )r — 2ay 2 = 0 , (14) 

and the resulting solutions for a are plotted as a function of r for various values of y in 
Fig.[9](a). For given a and y, the extrema occur at the r values where the line a = const, 
intersects the curve corresponding to a given y. Figure [9](b) shows the region of the 
phase diagram where the first-crossing probability is unimodal and where it is bimodal. 
Here the phase boundary is everywhere single valued, in contrast to the case of Fig. [4] 
the initial path is an unrestricted Brownian motion. As a result of the symmetry of 
F, the phase diagram is symmetric about the point x = k- Also, unlike in the case of 
unconstrained Brownian path, the phase boundary is single-valued. 




X X 

(a) (b) 

Figure 9. (a) Loci of the extrema of F{xx) bi the a-r plane for various values of 
X when the initial trajectory is a first-passage path to (X,T). For a given a and y, 
the extrema occur at the values of r where the line a = const, intersects the curve 
that corresponds to a given y. (b) Numerically obtained phase diagram from Eq. (lld|l . 
showing where in the parameter space the first-crossing probability, for a first-passage 
Brownian path from (0, 0) to ( X , T), is unimodal and where it is bimodal. The points 
marked a, b , c are at (0, § — 3\d2), (^, |), and (1, | — 3y/2), respectively. 
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3.2.2. Last-crossing probability We now determine when the particle crosses a specified 
intermediate level for the last time. In principle, we may perform this calculation by 
decomposing the trajectory into an initial segment from (0, 0) to the last crossing of 
x at time t, and a remaining first-passage segment from (x,t) to (X, T). The subtle 
feature here is that this second segment must also obey the constraint that this segment 
always remains greater than x, so that no additional crossings of x occur after time t. 
To satisfy the latter condition, we must also impose an absorbing boundary condition 
at x. Thus the particle begins the second segment at the absorbing boundary, so that 
the first-passage probability to (X, T ) would equal zero. To sidestep this pathology, we 
could start the particle at x+dx and take the limit dx at the end of the calculation. This 
limiting process is a delicate, however, and we therefore give an alternative approach 
that avoids any limiting processes. The price, however, is the necessity to break the 
trajectory into three segments (Fig. fTT)Tl . 



Figure 10. Zoom of the first-passage path in Fig. [I] with the decomposition of the 
segment from ti to T into a time-reversed first-passage path from (y,s) to (x,te) and 
a first-passage path from (y,s) to (X, T). 


In this alternative approach, the first segment is a Brownian path that starts at 
(0, 0) and reaches the last crossing of x at time tg. An absorbing boundary at X must 
be imposed to ensure that this segment never reaches X before T. We now break 
the remaining segment from (x,tf) to (X, T), into two sub-segments with respect to 
an arbitrary point (y,s), with x < Y < X and ti < s < T. The left sub-segment 
is a backward-propagating first-passage path from (y, s ) to (x, tf) and the right sub- 
segment is a forward-propagating first-passage path from (y,s) to (X, T). The choice 
of intermediate point (y, s ) is arbitrary, so that the final result must be independent of 
s after integrating over all y. 

Thus we have 


£(x, t ) 


Pa(x, t ) ff F A (y -x,s- t)F A (X -y,T- s)dy 
F(X,T) 


(15) 


Here F A (z,t ) denotes the first-passage probability in the presence of the absorbing 
boundary condition at X. Making use of standard results for the first-passage 
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Figure 11. (a) Comparison of the first-crossing probability F(x,t) and (b) the 

last-crossing probability £(y,r) for representative values of a for the specific case of 
X = 1/2. As a increases F(x, T ) remains symmetric, but changes from bimodal to 
unimodal, while £(x,t) merely moves systematically toward r = 1. 


probability in the interval [29], the integral in (fl5|l is 
[ F A (y-x,s-t)F A (X -y,T-s)dy 




nm sin 


m,n=1 


niry 

X-x 


sin 


m7l y ] e -^D[n\s-t)+mHT-s)\/(X-x) 2 dy 

X — x ) 


2tt 2 D 2 ^ 1W _ 

Tx^l^}- l)ne 

x 7 71=1 


n r ,2p-n 2 ir 2D {T-t)/( x - x ) 2 


(16) 


To obtain the third line, the orthogonality relation J 1 2 sin(n7r.x') sm(rmrx)dx = S mn has 
been used. 

Substituting the above expression, together with the previously derived forms for 
P A (x,t ) and F(X,T) in Eq. (fT5jh and using the dimensionless variables y, r and a, we 
finally obtain 


C(x, t ) = 


7 r 2 e Q 


v8a 2 - v /r(l —x) 3 


e ~ax 2 /r _ g— a(2 —x) 2 /t 


E(-d 


rc^g-nVp-rVIMl-x) 2 ] 


(17) 

When a -C 1, for 1 — r 3> ol (r not too close to 1), the first term in the sum dominates 
and the last-crossing probability reduces to 


e -A( 1-t) 

£{xx) -^ 17 — , (18) 

with A = 7r 2 /[4 ck(1 — x) 2 ]- I 11 Ibis limit, the last-crossing probability increases 

exponentially with r as r — y 1 (Fig. fill) . Conversely, when a 1, the last-crossing 
probability is sharply peaked at r = x- This feature is a result of the trajectory becoming 
nearly ballistic. When a 1, the first and last crossing times coincide because there is 
little stochasticity in the trajectory. 
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Thus far, we have tacitly assumed that the intermediate position is between 0 and X. 
With this constraint, every path from (0, 0) to ( X , T ) must necessarily pass through 
any intermediate position. We now investigate the first- and last-crossing probabilities 
when the intermediate position x < 0. ffere the probability that a first-passage path 
from (0, 0) to (X, T) actually reaches x is less than one. The decomposition of the first 
passage path into a segment from the start to the crossing point and from the crossing 
point to the final point still applies, but there are additional subtleties associated with 
the intermediate point being negative. 

Again, we separately consider the first- and last-crossing probabilities. It is 
convenient to define z — |x|, which is manifestly positive. The distribution of first- 
crossing times to z is formally given by 

„ F A (z,t)F(X + z,T-t) 

Z ’ } F(X, T) 

The constituent first-passage probability F^(z, t.) is subject to the constraint that X 
cannot be reached before time f; this requirement ensures that the full path is actually 
a first-passage trajectory to (X,T). The second segment from (z,t) to (X,T) is a 
first-passage path between these two points without any additional constraints. 

To obtain F^z^t), we again solve the diffusion equation in the interval between 
z = —x and X, with absorbing boundary conditions at both ends, and then compute 
the flux to the boundary point z for a particle that starts at the origin. By standard 
methods [29], this computation gives 



(19) 


„ , , 2ir D y. 

Fa{z ' 1) = 

Thus the first-crossing probability is 


n sm 


U7rZ 1 r -n 2 -K 2 Dt/(X+z) 2 

X + z 1 


( 20 ) 


2ttD 


T{z,t) = 


(X + z) 2 ^ 

v 7 n =1 


n sm 


n7TZ j -n 2 Tr 2 Dt/(X+z) 2 l-^ + ^l —(X+z) 2 /4D(T—t) 


X + Z 


V^TT D(T-t) 3 


\X\ 


3 -V 2 /4 DT 


V4nDT 3 

In dimensionless variables, the above expression simplifies to 


= 0 „ V1 , -3^j^nsm 


2a(l + x)(l-r) 3/2 ^ ' V! + X 


n7T * \ e -n 2 7T 2 r/[Aa(l+x) 2 } 


( 21 ) 


( 22 ) 


For a —> 0, for r>a, the first term in the sum is dominant so that 


~.—At 




(1 - r) 3 / 2 


Oi —y 0 , 


(23) 
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(a) (b) 

Figure 12. (a) The first-crossing probability P{x, T ) and (b) the last-crossing 

probability £(x,r) for y = —1/2 for a range of a values. The distributions have 
been rescaled so that they fit on the same range. 


with A = 7T 2 / [4ct(l + y) 2 ]. Conversely, for a —> oo, Fig. fl2l shows that J-{x, t) is peaked 
at r = The particle thus moves ballistically in the negative-x direction to reach 
—x in a fraction of the total time and then moves ballistically in the positive-x 
direction to reach X at time T. 

The distribution of last-crossing times has the identical form to that given in 
Eq. (11711 for the case where x > 0, except for the substitution x —> —x. For aCl, the 
last crossing necessarily occurs very close to the final time r = 1. For a > 1, the overall 
trajectory is a straight line in the negative-x direction to the intermediate point and 
then another straight line in the positive-x direction to the final point. In this limit, 
the first- and last-crossing probabilities become progressively more similar, as expected 
intuitively. 


5. Summary and Outlook 

We investigated the properties of intermediate crossings of Brownian and first-passage 
paths in one dimension that start at (0,0) and end at a final point (X,T). By using 
simple probabilistic arguments, we analytically determined the occupation probability 
at a general intermediate time, as well as the first-crossing and last-crossing times at 
a general intermediate location. These intermediate properties exhibit a number of 
surprising and anomalous features that depend, in an essential way, on a = X 2 /ADT. In 
the ballistic limit of a 1, a trajectory typically moves systematically from the starting 
point to the final point. In this case, one can infer intermediate-crossing properties by 
linearly interpolating the initial trajectory. 

The more interesting situation of the diffusive limit, where a<l, a first-passage 
Brownian path to ( X , T) must initially move away from this final point so as to not 
reach X before time T. The maximum excursion away from the starting point scales 
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as 1 jy/oL. The behavior of the first-crossing probability, the distribution of times 
when a path first crosses a specified intermediate point x is quite rich. When the 
ostensibly simpler example where initial path is a Brownian trajectory to (X, T) with no 
additional constraint, the first-crossing probability can be either unimodal or bimodal, 
with bimodality favored for small a and \ — x/X ss 0.75. Moreover, the phase boundary 
between unimodality and bimodality is, quite surprisingly, not single valued. When the 
initial path is a first-passage path to (X,T), the first-crossing probability again has a 
unimodal to bimodal transition and the corresponding phase diagram is symmetric in 
the a~x plane. 

Another unanticipated feature of the first- and last-crossing probabilities is that 
the results are invariant if there is an overall global drift in the Brownian trajectory. 
When one generalizes the expressions for the first- and last-crossing probabilities given 
in Eqs. (J4j), (ITTTh and (Tl5l) to incorporate a constant drift, all factors that involve the 
drift velocity cancel out. It would be worthwhile to understand the full ramifications of 
this simple observation. 

Finally, it is worth mentioning that the decomposition of the first-passage path 
in (fTTh provides the starting point for efficient simulations of intermediate crossing 
phenomena. The naive way to numerically determine intermediate crossing properties 
is by direct simulation of a random walk in which the times at which various 
intermediate positions are reached are recorded for a walk that reaches a given final 
point. From this data, one can reconstruct the first-crossing (as well as the last¬ 
crossing) probabilities. As a much more efficient alternative, one can first select a set of 
intermediate positions and then directly move a particle only between these intermediate 
positions. Correspondingly, the time between each of these macro-steps is incremented 
from the appropriate distribution of first passage times. This procedure can be made 
still more efficient by only allowing the walk to move to intermediate positions that are 
progressively closer to the final position. Thus arbitrarily long random walks can be 
simulated with a finite number of steps. 
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